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O ABSTRACT 



High resolution (0''15) Very Large Array observations of 7 mm continuum 
and H53a line emission toward the hypercompact H II region G28.20-0.04N re- 
veal the presence of large-scale ordered motions. We find a velocity gradient of 
10 3 km s -1 pc _1 along the minor axis of the continuum source. Lower resolution 
(l''0-2''3) radio recombination line observations indicate a systematic increase of 
line width from H30a to H92a. Under the assumption that the H30o; line does 
not suffer significant pressure broadening, we have deconvolved the contributions 
of statistical broadening (thermal, turbulent, and pressure) from large-scale mo- 
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tions. The pressure broadening of the H53a, H76a, and H92a lines implies an 
electron density of 6.9 xlO 6 cm -3 , 8.5 xlO 5 cm -3 , and 2.8 xlO 5 cm -3 , respectively. 
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^ ! 1. Introduction 

OO 

Hypercompact (HC) H II regions have about one tenth the size (<0.05 pc) and about 
a hundred times the density (n e >10 5 cm -3 ) of ultracompact (UC) H II regions. It has 
^ ■ been postu lated that H C H lis represent the transition from hot molecular cores to UC H 



II regions (IKurtzl 120001 ): that is, the transition from rapid accretion onto a massive central 



protostar to the termination of accretion and the emergence of a detectable radio H II region. 
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If HC H lis occupy this evolutionary niche, they should show evidence of ordered large-scale 
motions such as expansion, infall, bipolar outflows, rotation and shocks due to accretion 
and outflows. Also, owing to their high electron densities, recombination lines from the H + 
gas should suffer substantial pressure broadening in high quantum state radio recombination 
lines (RRLs). 

It is not clear if HC H II regions are simply smaller and denser versions of UC H II 
regions or if they are intrinsically different. Their relative sizes, however, suggest that stellar 
clusters may be responsible for UC H II regions, while HC H II regions may be produced by a 
single, or perhaps a binary, OB protostar. Some, but not all, HC H lis have unusually broad 
radio recombination lines (BRRLs) with full widths at half maximum intensity (FWHM) 
>40 km s" 1 . The typical FWHM of more evolved (UC and compact) H II regions is ~25-30 
km s _1 in all RRLs from principal quantum states n<100. The fact that not all HC H lis 
have BRRLs may imply that BRRLs are only present for a limited period during the earliest 
HC H II evolutionary phase when ordered large-scale motions are maximum. 



Sewilo et al.l (12004 and in prep.) observed the H92a, H76a, and H53a lines toward 



the HC H II region G28.20-0.04N (hereafter G28.20N) with resolutions of 2"3, 1'.'5, and l'/6, 
respectively, all of which are larger than the angular size of the HC H II region. They found 
FWHMs of 74 km s -1 , 58 km s _1 , and 40 km s _1 , respectively. The rapid change in line 
width with transition strongly suggests that pressure broadening is a dominant broadening 
mechanism in G28.20N. 

We have chosen G28.20N for a high resolution H53a study to determine the origin of the 
BRRLs in this source because: it is a bona fide HC H II region; it has already been observed in 
several RRLs at lower angular r esolution; several masers have bee n detected toward G28 .20N 
implying ext reme youth (H?Q - San et al.||l99~8l; OH 6035 MHz - Icaswell fc Vaile! ll995l: OH 



1667 MHz - lAreon et al.l l2000l : and, CH 3 OH - iMentenl Il99lh . The H53a line was chosen 



because it should be mostly free from pressure broadening, the nebula should be optically 
thin at this frequency, and the Very Large Array can achieve an angular resolution one tenth 
that of previous observations at 43 GHz. 
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2. Observations 

G28.20N was observed with the VLA0 at 7 mm (continuum and H53a line) with an 
angular resolution of 0''15. The observations were made on 2004 January 19 with the VLA 
in its B-configuration. At the H53a rest frequency of 42.95197 GHz, this provided an an- 
gular resolution of 0'.'17 x 0''13 at a position angle of — 4?8. A 25 MHz bandwidth was 
employed, with 31 spectral line channels of 781 kHz each. This provided a velocity reso- 
lution of 5.5 km s -1 , which was later Hanning smoothed to 11 km s -1 . Integration time 
on-source was about 1.4 hours. 3C 286 (1.47 Jy), 1851+005 (0.68±0.01 Jy), and 3C 345 
(4.81±0.18 Jy) were used as flux density, phase and bandpass calibrators, respectively. The 
data reduction, including self- calibration, followed the VLA high-frequency guidelines for 
calibration and imaging, using the NRAO package AIPS. The continuum image is shown in 
Figure 1. Continuum-subtracted line images were used to generate H53a maps. 



3. Results 

3.1. Morphology of the Source 

At our high angular resolution, G28.20N appears as a shell-like structure, extended 
in the SE-NW direction (Figure 1). There are two peaks, separated by 0''36 (2050 AU 



or 0.01 pc at a distance of 5.7 kpc; iFish et al.l 120031 ) . The brightness distributions along 



one-dimensional cuts through the source show a well-defined shell with low intensity at the 
center. The center of the shell-like nebula is at RA = 18 h 42 m 58?ll and Dec = -04°13'57"4 
(J2000). The source diameter at 10% of the peak (the 12a contour in Figure 1) is 9 S = 0"9 
(5100 AU; where 9 S = (9ra • #Dec) 1//2 )- The integrated flux density is 645±65 mJy, where the 
10% uncertainty reflects VLA calibration uncertainty in the 7 mm band. The synthesized 
beam brightness temperature using the peak brightness of 64 mJy beam -1 (located at the 
W side of the shell) is 1920 K, implying a peak optical depth of 0.3 based on T e = 7000 K, 
and a peak emission measure (EM) of 1.6 xlO 9 pc cm -6 . This EM implies an rms electron 
density of n e = 6.2 xlO 5 cm -3 . The synthesized beam brightness of the central cavity, 11 
mJy beam -1 , corresponds to 343 K, implying an optical depth of 0.05, an EM of 2.5 xlO 8 pc 
cm -6 , and an rms n P of 2.5 xlO 5 cm -3 . 



1 The Very Large Array (VLA) is operated by the National Radio Astronomy Observatory (NRAO), a 
facility of the National Science Foundation operated under cooperative agreement by Associated Universities, 
Inc. 
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3.2. Velocity Structure 

We examined the velocity structure of G28.20N by making contour plots of the H53a 
emission, using moment maps, position-velocity diagrams, and by Gaussian fitting of the 
line profiles at multiple positions. All these methods reveal a velocity gradient of 10 3 km s _1 
pc -1 , with velocities increasing from the NE to SW. 

In Figure 2 we show least-squares Gaussian fits to the H53a line integrated over different 
areas of the source, each approximately one synthesized beam in size. The positions are 
indicated in the continuum image of Figure 2. The boxes are located along the major and 
minor axes of the shell, with two boxes (B and D) sampling the regions of maximum emission. 
Box C covers the low intensity central region, while boxes F and I trace an extension in the 
SE-NW direction. We also analyzed the line emission integrated over the entire source (Box 
J). The parameters of the Gaussian fits are given in Table 1; Vlsr is the velocity relative 
to the local standard of rest, FWHMrf econt) is the line FWHM corrected for instrumental 
broadening, Sl is the flux density in excess of the continuum at line center, and Sc is the 
continuum flux density. In some directions (e.g. E and H) the line profiles are asymmetric 
with an additional velocity component apparent. All the asymmetries occur in the velocity 
range of ~1 10-130 km s _1 . Although weak, these features are real; we also detect a 130 km 
s _1 component in the H76a line. 

The sequence of five H53a line profiles shown in Figure 2 indicates a systematic shift in 
the central line velocity from the NE to SW. The dotted vertical line indicates the 88.2 km s _1 
systemic velocity of the nebula obtained from H53a data taken with 10 times lower resolution 
than shown in boxes A-E (upper right; HPBW ~ l."6)|E Thus, the centroid velocity of ~78 
km s^ 1 in Box A (NE side of G28.20N) is blueshifted and the centroid velocity of ~97 km s _1 
in Box E (SW side of G28.20N) is redshifted with respect to the systemic nebular velocity. 
The H53a lines suggest rotational motion with the axis of rotation in the SE-NW direction. 
There is no velocity pattern along the elongated SE-NW axis (i.e. no indication of an outflow 
along this axis). We suggest that the excess line broadening along this axis may be due to 
an outflow in the plane of the sky, expanding perpendicular to its axis of flow with a velocity 
of ~5 km s _1 . 



2 O ur 7 mm flux density and H53a line width (both B- and D-array) disagree with the results of lKeto et al 
(|2008h for their 7 mm/H53a data (S lnt =l.l Jy and FWHM = 33.4±1.5 km s" 1 ). To investigate this 
discrepancy, we re-reduced our lower resolution 7 mm data and the archival Keto et al. data (AK607) using 
consistent procedures for both data sets. We confirmed our results and found that the Keto et al. data 
are consistent with our own. For the AK607 data we find an integrated 7 mm continuum flux density of 
790±35 mJy and a FWHM for the H53a line of 38.7 ± 0.9 km s" 1 . 
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The velocity structure observed in G28.20N is consistent with a torus rotating with a 
velocity of 5.0 km s~ x at a radius of 0.005 pc (1030 AU; at the position of the peaks), which 
corresponds to a velocity gradient of about 10 3 km s _1 pc -1 or an interior mass of 28 M Q 
for Keplerian motion. 



3.3. Optical Depth Effects 



We compared integrated line profiles for all transitions and found that the H30a line 
velocity shifts by about 5-14 km s _1 with respect to the H92a, H76a, and H53a lines. The 
H92a, H76a, and H53a line velocity integr ated over the nebu la is 82.9±0.9 km s _1 , 78.9±0.8 
km s _1 , and 88.2±0.4 km s _1 , respectively. iKeto et al.l (120081 ) reports the H30a line velocity 
of 92.5±0.2 km s _1 . A similar trend of increasing RRL velocity with increasing frequ ency 



has been reported for other comp ac t H II regions, e .g., W 49A/ B (IDe Pree et al. 



W3(OH) (IBerulis & Ershovi Il983f ). IWelch fc Marrl (119871 ) and IKeto et al.l (119951 ) proposed 



19971) and 



that this is an optical depth effect, in which deeper layers with different velocities are probed 
at higher frequencies. Following their arguments, the optical depth decreases rapidly at high 
frequencies because the absorption coefficient is inversely proportional to the square of the 
frequency; thus the H92a, H76a, H53a, and H30a emission arise from different volumes 
of gas with H30a emission being more sensitive to the denser regions. They also argue 
that the velocity of the high frequency RRLs provides a good approximation to the average 
motion of the H II region provided that the random motions are the same for all densities 
throughout the nebula. The velocity of this innermost H II gas is expected to most closely 
indicate the stellar velocity; hence, for G28.20N, the H30a line velocity should provide the 
best estimate of the central star's radial velocity. The blueshifted velocities of the outer H 
II region, traced by the H92a, H76a, and H53a lines, imply expansion, consistent with a 
strongly over-pressured nebula. 



4. Line Broadening Mechanisms: Pressure Broadening vs. Dynamic Motions 



Keto et all (120081 1 detected the H30a line toward G28.20N using the SMA with an 



angular resolution of 1". The width of the H30a line is 20.9±0.6 km s _1 . At the high fre- 
quency of the H30a line, line width should be dominated by the thermal, turbulent, and 
large scale motions; pressure broadening is negligible due to the rapid decrease in pressure 
broadening with frequency. Thermal and turbulent broadening produce convolved Gaussian 
profiles whose center velocities generally coincide. Large-scale coherent motions can produce 
velocity-shifted profiles that may deviate from Gaussian shapes. The observed profiles de- 
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pend on the type and angular extent of the large-scale velocity components relative to the 
synthesized beam. In particular, as the angular resolution is increased, each resolution ele- 
ment should approach the combined thermal plus turbulent broadening of the gas, with its 
central velocity reflecting the dominant large-scale motion of the gas in that volume element. 
In some cases, two components of large-scale coherent motions may fall within a single beam, 
which may result in a double Gaussian. An example of this is an expanding bubble whose 
front and back faces are included in all observations or if the red and blue lobes of an outflow 
are projected onto the same resolution element. 

Molecular line observations of UC H II regions and massive cloud cores suggest that 
the typical molecular line width where HC H II regions form is ~5 km s _1 (FWHM; 



Cesaroni et al.l Il99ll . iPlume et al.l 119971 ). We adopt this value for the non-thermal (tur- 
bulent) contribution to the line width since the thermal width for molecules is <0.5 km s _1 . 
For 7000 K, the thermal contribution to the line width is 17.9 km s _1 (FWHM). Thus, the 
combined contribution (added in quadrature) of thermal and turbulent widths to a RRL 
profile is 18.6 km s _1 . We can determine the contribution from kinematic motions by sub- 
tracting this value in quadrature from the observed H30a line width. The residual after 
accounting for thermal and turbulent motion is 9.6 km s -1 . This excess represents the ap- 
proximate magnitude of the large-scale motions (detected in the H53a line) within G28.20N 
as the data is averaged over the entire source. For non-thermal widths of km s" 1 (no 
turbulent motions) and 10 km s _1 , the contribution from large-scale motions would be 10.8 
km s _1 and 4.1 km s _1 , respectively. The value of the non-thermal contribution to the line 
width is not essential to the following calculations of pressure broadening. 

The line width in km s _1 arising from thermal and turbulent motions is constant for all 
relevent transitions (from H30a to H92a). Thus, we can determine the influence of internal 
kinematics on the H53a, H76a, and H92a line widths. The subtraction in quadrature of 
the 9.6 km s" 1 kinematic broadening from the observed line widths (39.7±1.3, 57.6±2.2, 
and 74.4±2.6 km s" 1 ) gives FWHMs of 38.5, 56.8, and 73.8 km s" 1 for the H53a, H76a, 
and H92a lines, respectively; these line widths, "corrected" for kinematic broadening, still 
contain a pressure broadening contribution. 

To estimate the pressure broadening o f the H53a, H76a, and H92 a lines, we use an 



approximation for the Voigt profile width ([Gordon fc Sorochenkol 120031 . eq. 2.72), which 



relates the frequency widths of the Gaussian, Lorentz and Voigt profiles. The Gaussian width 
for low frequency lines is determined from the width of a high frequency line, unaffected by 
pressure broadening. Here we use the 16.167 MHz H30a line width. For the H53a line, the 
Gaussian width is 2.99 MHz and the observed line width is 5.69 MHz, implying a pressure 
broadening width of 4.03 MHz (or 28.1 km s _1 ). For the H76a line, the Gaussian width is 
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1.02 MHz and the observed line width is 2.82 MHz, implying a pressure broadening width 
of 2.43 MHz (or 49.5 km s" 1 ). For the H92a line, the Gaussian width is 0.58 MHz and the 
observed line width is 2.06 MHz, implying a pressure broadening width of 1.89 MHz (or 68.1 
km s _1 ). 

The pressure broadening width is proportional to n 7 4 T~ 01 n e ; thus, the pressure broad- 
ened lin es can be used to determine the true space electron density n e averaged over the 
nebula farocklehurst & Seatonl fl9~72l ) . Using T e of 7000 K, we find n e of (6.9±0.4) x 10 6 



cm" 3 , (8.5±0.4)xl0 5 cm" 3 and (2.8±0.1) x 10 5 cm" 3 based on the H53a, H76a, and H92a 
lines, respectively. The rms electron densities i nferred from the .7, 2, and 3.6 cm continuum 
observ ations are 1.7xl0 5 cm~ 3 , 9.3xl0 4 cm -3 (jSewilo et al.ll2004l . and in prep.), and 7.6xl0 4 



cm 3 (jSewilo et al.ll2004l ). respectively. The higher electron density from pressure broaden- 



ing is expected because it represents the actual space density of the gas that produces the 
broad lines as opposed to the rms density derived from the radio continuum. The fact that 
the average density decreases with higher quantum transitions also suggests that each line 
in the range 0.7 to 3.6 cm samples different volumes of gas, as also suggested by the change 
in line velocity with transition. 

The rms and true electron densities can be used to calculate the filling factor / for the 
nebula. Equating the Stromgren relation for a uniform density region (in terms of n rms ) with 
that of a clumpy region (in terms of n irne ) which is formed of N smaller clumps, one obtains 



Nr 



clump I ^"uniform 



true 



f. This is the ratio of the volume occupied by the dense 



gas to the total volume of the H II region; thus / is a measure of the degree of dumpiness 
of the nebula. Based on the 7 mm continuum and line data, we calculated the filling factor 
of 6xl0 -4 for the highest density gas we sample in G28.20N. 



5. Summary 

We have resolved the morphology and velocity structure of the HC H II region G28.20N. 
The continuum emission has a shell- like structure, with an inner radius of 1100 AU, an outer 
radius of 2500 AU, and an rms density contrast of 2-3 between the central cavity and the 
torus. 

Under the assumption that the H30a line is not significantly pressure broadened and 
that turbulence broadening is ~5 km s _1 , we have resolved the contributions of thermal, 
turbulent, and pressure broadening from large-scale ordered motions. We suggest that the 
ordered motion has two components: a velocity gradient of 10 3 km s" 1 pc^ 1 in the NE- 
SW direction along the minor axis which may indicate rotation of a torus around a 28 M 
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object; and line broadening (line splitting at some points) along the major axis which may 
indicate lateral expansion of an outflow along the major axis of the nebula. In this scenario, 
the outflow axis would have to be in the plane of the sky because line velocities along this 
axis are centered on the systemic velocity of the nebula. We suggest that the excess line 
broadening along the major axis is due to expansion of the outflow perpendicular to its axis. 

The line width from large scale motions averaged over the whole nebula is ~10 km s -1 . 
We find pressure broadening line widths for the H53o;, H76a, and H92a lines of ~28, ~50, 
and ~68 km s _1 , respectively. These spatially unresolved, pressure broadened lines were 
used to determine electron densities averaged over the nebula (not rms densities). The mean 
n e is 6.9 xlO 6 cm" 3 , 8.5 xlO 5 cm" 3 , and 2.8 xlO 5 cm" 3 based on the H53o;, H76a, and H92a 
lines, respectively. These densities are 4-40 times higher than the rms densities derived from 
continuum data. The fact that the average electron densities decrease with higher quantum 
transitions is consistent with the change in opacity with frequency. 

This study has for the first time resolved the contributions of thermal, turbulent, pres- 
sure, and ordered large-scale motions of BRRLs in a hypercompact H II region. The ob- 
servations show that large scale dynamics and high densities both contribute to BRRLs in 
G28.20N. This study was possible only because of the high angular resolution and sensitivity 
provided by the VLA. 

We thank Michael Burton, the referee, for useful suggestions that improved our manuscript. 
M. S. and E. B. C. acknowledge partial support from NSF grant AST-0303689. S. K. ac- 
knowledges partial support from DGAPA, UNAM grant IN106107. 
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Table 1. Gaussian Fits to the High Resolution H53a RRLs a 



Box 


Vlsr 


FWHM dcconv 








(km s _1 ) 


(km s- 1 ) 


(mJy) 


(mJy) 


A 


77.6 ± 1.7 


34.9 ± 4.9 


5.1 ± 1.0 


12.1 ± 0.6 


B 


85.8 ± 0.5 


33.8 ± 1.3 


22.5 ± 1.1 


53.4 ± 0.6 


C 


92.1 ± 0.7 


35.9 ± 2.2 


6.7 ± 0.6 


23.2 ± 0.6 


D 


95.8 ± 0.4 


31.6 ± 1.0 


26.4 ± 1.1 


66.5 ± 0.6 


E b 


r 97.3 ± 1.0 


28.9 ± 3.0 


7.8 ± 1.0 


17.9 ± 0.6 




1127.9 ± 3.0 


9.4 ± 5.0 


1.8 ± 1.3 




F 


89.9 ± 1.3 


29.6 ± 3.7 


4.0 ± 0.7 


10.4 ± 0.6 


G 


87.7 ± 0.8 


36.6 ± 2.4 


11.4 ± 1.0 


29.5 ± 0.6 


H b 


r 82.8 ± 1.7 


26.5 ± 3.3 


10.1 ± 1.7 


34.7 ± 0.6 




1110.3 ± 4.0 


21.1 ± 7.4 


3.6 ± 1.8 




I 


83.2 ± 2.1 


28.1 ± 5.9 


3.4 ± 0.9 


8.5 ± 0.6 


J 


89.5 ± 0.6 


39.8 ± 1.7 


238 ± 14 


641 ± 4 



a Spectra were Hanning smoothed to a resolution of 11.0 km 
s -1 . Uncertainties are the formal la (68.3% confidence level) 
deviations of the data from the Gaussian fits. 

b Two Gaussian components were required to adequately fit 
the observed line profile. 
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G28.20-0.04 N 7 mm 




18 42 58.20 58.16 58.12 58.08 58.04 

RIGHT ASCENSION (J2000) 



Fig. 1. — The 7 mm image of G28.20— 0.04 N. The continuum intensity peak is 64.0 
mJy beam" 1 . The contours are -3, 3, 6, 9, 12, 20, 30, 50, 70, 90, 110 x 0.53 mJy beam^ 1 
(the image rms). The synthesized beam of 0''17 x0''13 is shown in the lower left corner. 
The insets show intensity cuts through G28.20N; position angles indicate rotation East from 
North. 




Fig. 2. — A sequence of five H53o; line profiles shown in the diagonal direction illustrates 
a shift in the line central velocity of ~20 km s" 1 along the NE-SW direction. The spectra 
were Hanning smoothed to a resolution of 11.0 km s _1 . The dotted vertical line indicates the 
systemic velocity of the H53a line profile integrated over the entire nebula from data with 10 
times lower resolution (upper right; HPBW ~ 1"6, VLA D-Array; S int = 711 mJy). The lines 
were integrated over rectangular areas indicated in the 7 mm continuum image (lower left). 
The corresponding box letters are given at the top of each panel. In the lower left corner of 
the continuum image, the size of the integration box is compared with the synthesized beam 
of 0''17 x 0''13; the two have very nearly the same area. 



